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The measurements of the Higgs boson and top quark masses can be used to ex- 
trapolate the Standard Model Higgs potential at energies up to the Planck scale. 
Adopting a NNLO renormalization procedure, we: i) find that electroweak vacuum 
stability is at present allowed, discuss the associated theoretical and experimental 
errors and the prospects for its future tests; ii) determine the boundary conditions 
allowing for the existence of a shallow false minimum slightly below the Planck scale, 
which is a stable configuration that might have been relevant for primordial infla- 
tion; iii) derive a conservative upper bound on type I seesaw right-handed neutrino 
masses, following from the requirement of electroweak vacuum stability. 



I. INTRODUCTION 

The recent discovery of a particle consistent with the Standard Model (SM) Higgs 
boson, announced by the ATLAS [1] and CMS [2] collaborations at CERN, is a milestone 
in particle physics; adding in quadrature statistical and systematic errors, the mass of the 
particle turns out to be in the range 124.8 — 126.5 GeV at 2a. 

Here we assume that the new particle is actually the SM Higgs boson and study the 
implications that its mass value, together with other relevant parameters such as the top 
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quark mass and the strong gauge coupling, have on the behavior of the Higgs potential at 
very high energy scales and, in particular, for the sake of electroweak vacuum stability. 

The project of extrapolating the Higgs potential up to the Planck scale is a long 
standing one [3-5], and was revamped in the fall of 2011 [6-9] after the first LHC hints 
of a Higgs boson [10]. Recently, the tools for a Next-to-Next-to-Leading Order (NNLO) 
renormalization procedure were derived [11—14]. So, there are now all the ingredients to 
carry on this long standing project. Clearly, the extrapolation is based on the assumption 
that there is a desert up to the Planck scale or, better, that possible new physics do 
not affect significantly the running of the Higgs quartic coupling, dominating the Higgs 
potential at high energy. 

It is interesting that the recently discovered experimental Higgs mass range, combined 
with the experimental top mass range, indicates a particularly intriguing high energy 
behavior of the Higgs potential, close to the transition between electroweak vacuum sta- 
bility and metastability. This is due to the fact that, for these Higgs and top mass values, 
the Higgs quartic coupling can be very small or even negative. Since the dependence on 
the top mass is strong and quite subtle, it is not surprising that different groups slightly 
disagree in the interpretation of the results, some of them favoring [13] and some other 
disfavoring [14] electroweak vacuum stability. 

Traditionally the top pole mass was used in the analysis, but it has been pointed out 
[15] that the top pole mass value used in previous analyses and taken to be the one 
measured at Tevatron, mf ? = 173.2 ± 0.9 GeV [16], is not unambiguously derived and 
that a more careful derivation should be rather based on the running top mass in the MS 
scheme, ml{rrit) = 163.3±2.7 GeV. As shown in [15], the top pole mass range consistently 
derived from the running one, m t = 173.3 ± 2.8 GeV, is plagued by a larger error than 
the Tevatron measurement considered in [14], rescuing electroweak vacuum stability. 

In our analysis we keep as a free parameter the running top mass, rather than the pole 
one. In this way we completely avoid the theoretical uncertainties associated to the top 
Yukawa matching procedure. As we are going to discuss, the theoretical error associated 
to the Higgs quartic coupling matching [13, 14] turns out to be smaller than the one 
induced by the experimental uncertainty on the strong gauge coupling, a 3 (mz)- Given 
the above mentioned range for the running top mass [15], we find that electroweak stability 
is allowed in the whole Higgs mass range [1, 2]. Stability could be soon excluded if values 
of the running top mass rnt(m t ) < 163 GeV could be excluded by the LHC. Otherwise, 
testing electroweak vacuum stability would become very challenging, since this would 
require precision measurements of the Higgs and top masses, and also of a^{mz). 

A stable Higgs potential configuration which deserves particular interest is a shallow 
false minimum close to the Planck scale, which could have been relevant for primordial 
inflation [7, 17, 18]. We show that such configuration is realized only if the Higgs quartic 
coupling and its derivative satisfy very specific boundary conditions, possibly having a 
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deep origin in quantum gravity. 

As well know, new physics in addition to the SM is required to explain neutrino masses 
and mixings, and also dark matter. The mechanism responsible for neutrino masses could 
affect the Higgs quartic coupling; as an example, we consider the impact that the inclusion 
of neutrino masses via a type I seesaw has on electroweak stability, discussing in some 
detail the shallow false minimum configuration. 

The paper is organized as follows. In sec. II we discuss the input parameters and 
the NNLO renormalization procedure used to extrapolate the Higgs potential up to the 
Planck scale. An analysis of electroweak vacuum stability and the associated constraints 
on the top and Higgs masses, with a detailed discussion of the theoretical errors and the 
prospects for the future, are presented in sec. III. In sec. IV we investigate the boundary 
conditions leading to the particularly interesting configuration of a shallow false minimum 
below the Planck scale. Sec. V is devoted to the upper bound on the seesaw right-handed 
neutrino masses following from the requirement of electroweak vacuum stability. Con- 
clusions are drawn in sec. VI. Appendix A contains the relevant formulas for the NNLO 
running procedure in the SM and, in appendix B, those to incorporate the type I seesaw 
mechanism. 

II. INPUT PARAMETERS AND RENORMALIZATION AT NNLO 

The normalization of the Higgs quartic coupling A is chosen in this paper so that the 
potential for the physical Higgs 4>h contained in the Higgs doublet H = (0, (0# + v)/y/2) 
is given, at tree level, by 



where v = l/(v / 2G M ) 1/2 = 246.221 GeV and = 1.1663787(6) x KT 5 /GeV 2 is the Fermi 
constant from muon decay [19]. The approximation in eq. (1) holds when considering large 
field values. According to our normalization, the physical Higgs mass satisfies the tree 
level relation m 2 H = Xv 2 /3. In addition, the mass of the fermion / reads, at tree level, 



In order to extrapolate the behavior of the Higgs potential at very high energies, we 
adopt the MS scheme and consider the Renormalization Group (RG) evolution for the 
relevant couplings which, in addition to the Higgs quartic coupling A, are the gauge g, 
g', g 3 , and the top Yukawa h t couplings. We work at NNLO, namely 3-loops for the 
/3-functions and 2-loops for the matching conditions at some suitable scale. 

It is customary to introduce the dimensionless parameter t = log [i/rnz, where \x stands 
for the renormalization scale and mz is the Z boson mass. The RG equations for the 
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relevant couplings are then given by 

d 



dt X(t) = K$ + #fif> + K 3 /?f 
j t 9(t) = kP® + Jfif + J {if 

f t 9'(t) = K$ + K^f + .^f 



(2) 



where k = l/(167r 2 ) and the apex on the /3-functions represents the loop order. The 1- 
loop and 2-loop expressions for the /3-functions can be found e.g. in ref. [20] (see also [21- 
26]). Recently, the complete 3- loop /3-functions for all the SM gauge couplings have been 
presented by Mihaila, Salomon and Steinhauser in ref. [11], while the leading 3- loop terms 
in the RG evolution of A, h t and the Higgs anomalous dimension have been computed 
by Chetyrkin and Zoller in ref. [12]. For the sake of completeness, the expressions for the 
/3-functions up to 3-loops are collected in appendix A. 

The matching of the running gauge couplings is done at the Z boson pole mass, mz- 
The numerical values used for the related MS observables are taken from the latest Particle 
Data Group SM fit results [19]: 

a -^(m z ) = 127.944 ± 0.014 , a 3 (m z ) = 0.1196 ± 0.0017 , (3) 
sin 2 e w (m z ) = 0.23116 ± 0.00012 , m z = 91.1874 ± 0.0021 GeV . 

To match the MS running quartic coupling A(/i) with the Higgs pole mass m# is more 
complicated and requires to exploit an expansion, 

*(/*) = E A(n) M = 3 ^ { 1 + 5 hM+ 4 2) 0*) + -) > (4) 

n=l,2,3,... 

which is know at present at NLO: 8^ (/i) is the 1-loop 0(a) result of Sirlin and Zuc- 
chini [27] while S H (//) is the recently calculated 2-loop result, composed by a QCD con- 
tribution of 0(aa 3 ) [13, 14] and a Yukawa contribution [14]. More details can be found in 
appendix A. As well known, there is some arbitrariness in the choice of the matching scale 
H in eq. (4), which introduces a "theoretical" error in the RG procedure. In this work, we 
choose to perform the matching of the Higgs quartic coupling A at the scale \i = m#. The 
theoretical uncertainty is estimated by performing the matching also at different scales 
and by evolving A via RG running until \x = ma- The spread in the numerical values 
obtained for A(m#) can then be used to infer the magnitude of the theoretical error. 

This is illustrated in fig. 1, assuming for definiteness a top pole mass m t = 172 GeV. The 
dashed and solid curves show the value of A(m#) obtained by including the corrections up 
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to 1-loop and 2-loop respectively, for various choices of the matching scale: from top to 
bottom fM = mz, ma, rnt, 2m//. One can see that, working at the 1-loop, the theoretical 
uncertainty is about 5%. The inclusion of the 2-loop corrections given in ref. [14] reduces 
the theoretical uncertainty down to about 0.7%. Notice also that the preferred region 
shrinks to small A values and that /i = mz and fi = m# nearly overlap. More generally, 
one can use the following expression for the 2-loop result, 

\{m H ) = 0.8065 + 0.0109 (m H [GeV] - 126) + 0.0015 (m t [GeV] - 172) ±°;°°° 2 , (5) 

where the mean value refers to fi = mn- 



0.83 - 
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FIG. 1: Value of A(m#) obtained by performing the matching at different scales fi, indicated by 
the labels, as a function of mn- The solid (dashed) lines are obtained by including corrections 
up to 2-loop (1-loop). We fixed rri( = 172 GeV (for different values see eq.(5)). 

It is common to extrapolate the MS top Yukawa coupling h t (/j,) from the matching 
condition between the running top mass m7(/z) and the top pole mass m t : 

thO*)A= = rfh(p) = m t {l + StQij) , W = 5f M + 5? ED (ri + 5? CD ^) , (6) 

where + represent the electroweak contribution, which is known at 1-loop [28] , 
while 5? CD is the QCD one. The QCD 1-loop result is known since many years [28]; 
the QCD 2-loop and 3-loop results as a function of the matching scale \x are given in 
[29] (see also [30-34]). The matching is usually done at the top pole mass scale, and the 
theoretical error associated to the arbitrariness of the matching scale can be estimated 
as before, namely by comparing the values of h t (m t ) obtained with different matching 
scales. This is represented in fig. 2, where the curves are obtained by working at 2-loop 
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and using, from bottom to top, \i = m Z} m t ,2m t . The plot shows that the associated 
theoretical uncertainty is about 2%. The analytical expression for h t {m t ) is: 

h t {m t ) = 0.933 + 0.006 (m t [GeV] - 172) t ll . (7) 

The procedure adopted in previous analyses of the stability of the electroweak vacuum 
was to use the experimental value of mt, identified with the one measured at the Tevatron 
by the CDF and DO collaborations, m e t xp = 173.2±0.9 GeV [16], to extrapolate the running 
Yukawa h t (m t ) via eq. (7). However, as discussed in ref. [15], it is not meaningful to use 
the mass parameter provided by the Tevatron as the pole top mass to be inserted in 
eq. (7): the running top mass in the MS scheme is instead a well defined parameter that 
can be directly extracted at NNLO from Tevatron measurements of the inclusive top pair 
production cross-section, giving rnl{m t ) = 163.3 ± 2.7 GeV [15]. So, it is conceptually 
more robust and practically more convenient to extract the top Yukawa coupling directly 
from fnl(m t ), as will be done in the following 1 . Our results will thus be presented as a 
function of m^(mj). 

Notice that, according to eq. (7), the value of the top pole mass can be easily recov- 
ered via the relation m t = mr t {m t ) + 9.6 jJJ GeV, which however is plagued by a large 




FIG. 2: Values of ht{mt) and fnt(mt) as a function of mj. The curves are obtained by matching 
at different scales, which are indicated by the labels. We fixed mjj = 126 GeV for definiteness 
but the results do not significantly dependent on run, provided it is chosen in its experimental 
range. 



At difference, ref. [15] proceeds in a more complicated way: the value of fn^(m t ) is translated into a 
value of m t , to be inserted in the expression of the lower bound on m# ensuring electroweak vacuum 
stability as derived in ref. [14]. 
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uncertainty. In ref. [15] it was found that, by doing a scheme transformation to NNLO 
accuracy from the running to the pole top mass, the range m^(m t ) = 163.3 ± 2.7 GeV is 
equivalent to mt = 173.3 ± 2.8 GeV. Hence, while displaying our results as a function of 
fn t {m t ) as already stated, motivated by the results of ref. [15], in some plots (as the one 
in fig. 5) we will link the value of the top pole mass to the running mass via the simple 
relation m t = WT t (m t ) + 10 GeV. 

Before presenting the results of our analysis in the following sections, we recall that, 
in order to carefully study the shape of the Higgs potential at high energy, one should 
consider the renormalization improved effective potential. This can be done by introducing 
an effective coupling, A e //(/i) = A(/i) + AA(/i), so that 

V.fMn) = ■ (8) 

The expression for AA(/i) is known up to 2-loop [4, 20] (and given, for instance, in [14]). 
Since the scalar contribution is not well defined when A is negative, in the following we 
consider the renormalization improved potential at the tree level, and identify fi with <pH- 
It is well known that this simplification has a negligible impact in the determination of 
the vacuum stability bound to be discussed in the next section. 



III. ELECTROWEAK VACUUM STABILITY 

The experimental region of the values of the Higgs and top masses is very intriguing 
form the theoretical point of view, since the Higgs quartic coupling could be rather small, 
vanish or even turn negative at a scale slightly smaller than the Planck scale. Accord- 
ingly, the behavior of the Higgs potential at high energy changes drastically: if A(^) is 
always positive, the electroweak vacuum is a global minimum, possibly accompanied by 
another local minimum just below the Planck scale, which could have played a role in 
primordial inflation [7, 17, 18]; if A(/x) turns negative below M Ph the electroweak vacuum 
correspondingly becomes metastable [4, 5]. 

These drastically different possibilities for the behavior of the renormalization improved 
Higgs potential at high energy are illustrated in the left plot fig. 3, where m# = 126 GeV 
and some specific values for mt(mt) have been selected, increasing from top to bottom. 
The right plot shows the associated values of A(/z). Let start considering the value 
rn t {m t ) = 161.989 GeV. Increasing the latter by just 1 MeV, the potential develops an 
inflection point; notice that the associated A(/z) becomes as small as O(10~ 5 ). Increasing 
again fnt(m t ) by about 200 keV, the minimum of A(/i) is equal to zero: a second vacuum 
degenerate with the electroweak one is obtained. Further increasing rfr t (mt) makes A(/i) 
turn negative: the electroweak vacuum becomes metastable. 

The dashed curve in the right plot in fig. 3 shows the evolution of j3\ (/x) = d\(fi)/dt 
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H [GeV] // LGeVJ 

FIG. 3: The SM Higgs potential (left) and the quartic Higgs coupling (right) as functions of the 
renormalization scale [/,, for ran = 126 GeV and different values of rni(mt), increasing from top 
to bottom by the amount indicated by the labels. The dashed curve in the right plot shows the 
associated value of (3\(n). The other input parameters are fixed at the central values discussed 
in the previous section. 

for the same parameter values; there is only a single dashed curve because Px(fi) mildly 
depends on m~ t {m t ) if the latter is in the range 161 — 163 GeV. Let call /ig the renor- 
malization scale such that /3\([ip) = 0. Clearly, only in the case of two degenerate vacua 
the conditions Px^Hp) = and \(fj,p) = are simultaneously met. For a shallow false 
minimum we instead have /3\(/j,p) = and \{np) = O(10~ 5 ), as already mentioned. 

In fig. 4 we show how \x$ depends on rnt(mt), for various values of run- It is remarkable 
that fip is maximized and nearly constant for the values of fni(m t ) for which A(/i) is very 
small. 



6x io 17 ; 
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FIG. 4: The scale as a function of mt(mt) and for different values of m#, as indicated by 
the labels. 
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We now turn to the determination of the points in the plane [m#,m^(m t )] allowing for 
the existence of a second minimum degenerate with the electroweak one. These points 
belong to a line separating the stability from the metastability region, see fig. 5: in the 
lower part of the plot A(/x) is always positive, while in the upper part it becomes negative 
before reaching the Planck scale. The configuration of a shallow false minimum belongs 
to the stability region, but the associated points are so close to the transition line that 
they could not be distinguished by eye. 
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FIG. 5: The solid (black) line marks the points in the plane [run , nTt(mt)] where a second 
vacuum, degenerate with the electroweak one, is obtained just below the Planck scale. The 
(red) diagonal arrow shows the effect of varying a^{mz) = 0.1196±0.0017; the (blue) horizontal 
one shows the effect of varying (the matching scale of A) from mz up to 2m h. The shaded 
(yellow) vertical region is the 2a ATLAS [1] and CMS [2] combined range, ran = 125.65 ± 0.85 
GeV; the shaded (green) horizontal region is the range rfTt{mt) = 163.3 ± 2.7 GeV, equivalent to 
m t = 173.3 ±2.8 GeV [15]. 

The transition line of fig. 5 was obtained with the input parameter values discussed in 
the previous section and by matching the running Higgs quartic coupling at ran- Clearly, 
it is also important to estimate the theoretical error associated to experimental ranges of 
the input parameters and the one associated to the matching procedure. To illustrate this, 
we consider in particular the point on the transition line associated to the value m# = 126 
GeV; for such point, A and (3\ both vanish at a certain scale [ip (see fig. 4). The arrows 
show how, if some inputs or the matching scale are changed, the position of this point 
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have to change in order to keep having, at the same scale [ip, a vacuum degenerate with 
the electroweak one. The diagonal arrow is obtained by varying the strong coupling in 
its allowed range, (mz) = 0.1196 ± 0.0017; the short (long) dashed line shows how the 
solid line would move if a^ivrtz) were equal to its minimum (maximum) presently allowed 
value. Notice that the error on a^(mz) induces an uncertainty in both the Higgs and top 
masses of about ±0.7 GeV. The variation of the other input parameters induces a much 
smaller effect. The horizontal arrow is instead obtained by varying /ia, the matching scale 
of the Higgs quartic coupling, from fi = m z to fi = 2m H ; notice that the associated error 
is very asymmetric (see fig. 1): essentially it can only enhance m#, by at most 0.5 GeV. 
Clearly, similar considerations apply to each point of the transition line. 

Since stability can be achieved in the whole experimental range for m# (shaded vertical 
region), but this is not the case for m~l(m t ) (shaded horizontal region), it is more convenient 
to write down the condition of electroweak vacuum stability under the form of an upper 
bound on the top mass: 

m£(m t )[GeV] < 162.0 + 0.47 (m^[GeV] - 126) ± 0.7 aa - 0.2 MA . (9) 

The first error is associated to the experimental error on a^[m z \ while the second ac- 
counts for the theoretical error induced by the matching of A. Since the latter is smaller 
and can be easily accounted for by just shifting mu upwards by at most 0.5 GeV, it will 
not be explicitly shown in the following plots. 

Given the present experimental situation, in order to understand whether we live in 
a stable or metastable vacuum it is crucial to better determine m~t(m t ). As discussed in 
[15], after LHC the Higgs mass will presumably be known with an accuracy of (9(100) 
MeV [35], but the precision on the top mass would improve only by a factor of two. 
For instance, if the whole range of m~t(mt) below 163 GeV would be excluded, we would 
conclude that our vacuum is metastable; otherwise the investigations should continue. 

A self-consistent and precise determination of the top quark mass can best be per- 
formed at a high-energy electron-positron collider, with a planned accuracy of (9(100) 
MeV. Moreover, at an electron-positron collider a^mz) could be determined with an 
accuracy close to or better than Aa^(m z ) = 0.0007 (this precision is sometimes currently 
adopted [13, 14] but cannot be considered to be conservative according to ref. [15]). At 
this stage, if the stability region will still have an overlap with the allowed ranges of the 
top and Higgs masses, we will be mostly limited by the theoretical uncertainty associated 
to 

Notice also that it is not realistic to hope to distinguish the case of two degenerate 
minima with the one of a shallow false minimum, since the difference in the top mass is 
just about 200 keV (see fig. 3). 
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IV. SHALLOW FALSE MINIMUM 

It is interesting to study in some detail the boundary conditions which must be satisfied 
in order to have a very shallow false minimum just below the Planck scale, since it could 
be relevant for inflation [7, 17, 18]. 

To study this particular configuration, we denote with \l% the renormalization scale 
where the Higgs potential has an inflection point; we also recall that \ip has been defined 
to be the scale where A(/i / g) = and /3\(jip) = are simultaneously fulfilled. Both fa 
and \xp increase 2 with m#, as shown in fig. 6, where the shaded region accounts for the 
experimental range of otz{mz)- 
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FIG. 6: Values of [a (left) and fj,p (right) as a function of ran- For the solid lines, the input 
parameters are fixed at their central values and the matching of A is done at fi = ran ■ The shaded 
regions shows the uncertainty induced by the experimental error of a^{mz) = 0.1196 ± 0.0017: 
the short and long dashed curves refer to the lower and upper value at la, respectively. 

It is interesting that, for the whole experimental range of m#, a shallow false minimum 
is obtained only if the following boundary condition holds: 

\{Hp) ~ (8.75 ± 0.15) x 10~ 6 . (10) 

One could speculate that such a special value could be dictated by an eventual merging 
with quantum gravity [6, 36]. In the left plot of fig. 7 we show that nas a mild 

dependence on in the right plot we show instead the value of the Higgs potential at 
the inflection point, which turns out to be of (9(10 16 ) GeV. As before, the shaded regions 
account for the experimental range of ^(m^). 

As pointed out in [17], a way of testing the hypothesis that inflation occurred when 
the Higgs field was trapped into a shallow false vacuum below the Planck scale is to look 



2 Notice that /ij is slightly smaller than fig. This can be easily understood, since the condition for having 
an inflection point at ^ reads (3\(fj,i) = —AX(fj,i) < 0, which implies /ij < \ip. 
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FIG. 7: Left: the value of A(/xg) as a function of m#. Right: the Higgs potential at \i{ and the 
associated prediction for r as a function of itlh- The short and long dashed curves refer to the 
la lower and upper values of a^{mz), respectively. The shaded (yellow) vertical region marks 
the preferred range of m# at 2a [1, 2]. The upper region in the right plot is excluded because 

< 



r^0.2 [37]. 



at the tensor-to-scalar ratio r of cosmological perturbations. The amplitude of density 
fluctuations in the observed Universe as seen by the CMB and Large-Scale structure data 
is parametrized by the power spectrum in /c-space, P s (k) = A 2 R (k/k ) ns ~ 1 , where is 
the amplitude at some pivot point ko, whose best-fit value is A R = (2.43±0.11) x 10~ 9 at 
k = 0.002 Mpc -1 [37]. In models where inflation happened while the Higgs was trapped 
in the shallow minimum [7, 18], the Higgs potential at the inflection point and the amount 
of gravity waves that can be produced - parametrized via the tensor-to-scalar ratio r - are 
linked via a simple relation: 

2 IV (ja) 
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37r 2 r M 4 ' 

where M is the reduced Planck scale. Such prediction for r is reported in the right plot 
of fig. 7. Notice that, for these models, only if m# is in its upper allowed range and 
otz{ m z) is quite low, there are chances for the Planck satellite mission [38] to measure r. 
However, the forthcoming experiment EPIC [39] should be able to test r down to 10 -2 , 
while COrE [40] down to about 10~ 3 . 



V. CONSTRAINTS ON THE SEESAW MECHANISM 

We now consider the effect of including neutrino masses via a type I seesaw. This issue 
has been already considered in a series of papers [8, 41-44]. 

Although the precise amount of the effect is quite model dependent, here we obtain 
a conservative estimate of the effect by considering only one right handed neutrino with 
mass M u , associated to a light Majorana neutrino with mass m v = 0.06 eV, the scale of 



13 



the atmospheric oscillations. This is supported by the following argument. 

It is well known that the /3-function of the Higgs quartic coupling is affected only 
if h u (n), the Yukawa coupling of the Dirac mass term (defined only for fi > M u ), is 
large enough. As the top Yukawa coupling, also the neutrino Yukawa coupling induces a 
suppression of the Higgs quartic coupling at high energy. By increasing M v and m u , the 
neutrino Yukawa coupling at the threshold scale M v also increases: 

M M„H2 V ^ffi. (12) 

This justifies that the fact that we equate ra v to the the atmospheric mass scale, about 
0.06 eV, which is the lowest possible value for the heaviest among the three light neu- 
trinos. In addition, two other Majorana neutrinos with masses lighter than m v can be 
accommodated via the seesaw but, if their right-handed neutrinos are lighter than M v , 
the associated Dirac Yukawa couplings are naturally expected to be smaller, and their 
effect on A(/i) negligible. 

In Appendix B we provide the additional terms (with respect to the pure SM) for the 
relevant /3-functions, above and below the scale M v . 

Since the effect of h v is a suppression of A, a configuration with a stable electroweak 
vacuum in the SM, could be rendered metastable because of the addition of the seesaw 
interaction. For a fixed value of m#, one can find the upper bound on M v following 
from the requirement that the electroweak vacuum is not destabilized. As shown in fig. 8 
for nin — 126 GeV (but similar upper bounds are obtained in the whole experimental 




FIG. 8: Upper bound on M u as a function of the running top mass, following from the require- 
ment that the electroweak vacuum is not destabilized because of the inclusion of the seesaw, for 
mu = 126 GeV. The shaded region is obtained by varying 03(771^) in its la range. 
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range of m H ), such upper bound strongly depends on the top mass 3 and is affected by an 
uncertainty which is mainly due to oi^imz) (shaded region). The smaller the top mass 
is, the more the configuration is stable and the less stringent is the M v upper bound, 
M v < 3 x 10 14 GeV. But increasing the top mass, the electroweak vacuum becomes less 
stable and the upper bound on M v becomes accordingly more and more stringent. 

Let consider in particular the upper bound on M v needed to avoid destabilization of an 
inflection point configuration, as the one depicted via the dashed line in fig. 9. Notice that 
an inflection point becomes a not so shallow local second minimum if M v ~ 10 11 GeV and 
that electroweak vacuum destabilization is avoided only if the condition M v < 2 x 10 11 
GeV is satisfied. The latter bound might be relevant for models of inflation based on the 
SM shallow false minimum [7, 17, 18]; note however that it is well compatible with the 
thermal leptogenesis mechanism to explain matter-antimatter asymmetry, for which the 
lower bound on the lightest Majorana neutrino is about 5 x 10 8 GeV [45]. 

Clearly, the neutrino Yukawa coupling y v is not the only additional term beyond the 
SM capable of modifying the running of A at high energy. Always in the context of type I 
seesaw, in the case that the vacuum expectation value of a singlet scalar field S (violating 
the lepton number by two units) is actually at the origin of the right-handed Majorana 
neutrino mass, the 5* couplings induce an enhancement of A, thus helping the stability of 




[GeV] 

FIG. 9: The Higgs potential as function of the renormalization scale, for mn = 126 GeV and a 
value of the top mass leading to an inflection point configuration in the SM case (dashed curve) . 
The lower curves display the effect of adding the seesaw, with three increasing values of M v 
from top to bottom (solid curves). 



3 This dependence was not considered in the previous literature. 
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the electroweak vacuum [46]. Such effect is indeed generically expected when adding to 
the SM a singlet field S [46, 47]. 

VI. CONCLUSIONS 

The recent discovery of a particle consistent with the SM Higgs boson [1, 2] provides 
a strong motivation to pursue [6-9, 13-15, 17] the old project [3-5] of investigating the 
behavior of the Higgs potential at very high energies. In particular, one would understand 
whether the electroweak vacuum is a global minimum up to the Planck scale, namely 
whether we live in a stable vacuum. In particular, a stable configuration which deserves 
a special interest is a shallow false minimum below the Planck scale: the Higgs field 
could have been primordially trapped there, leading to a stage of inflation [7, 17, 18]. 
Stability below the Planck scale is required also in Higgs inflation models with non minimal 
gravitational couplings [13, 48]. 

In our analysis, we adopted the recently derived tools for a NNLO renormalization 
procedure [11-14]. At difference of previous analyses, we considered as free parameter the 
running top mass rather than the (Tevatron) top pole mass, as suggested in ref. [15]. 

Given the present range of the running top mass and of the Higgs mass, we found 
that electroweak vacuum stability is at present allowed, as shown in fig. 5. To further 
test stability, a more precise measurement of the top mass would be crucial. As apparent 
from the stability condition of eq. (9), in case that LHC will not exclude values of the 
running top mass below 163 GeV, an electron-positron collider would probably be needed 
to discriminate between stability and metastability. 

We also determined the high scale boundary conditions allowing for a shallow false 
minimum slightly below the Planck scale, A(/i ( g) ~ 10~ 5 (fip is the renormalization scale 
were the /3-function of the Higgs quartic coupling vanishes), and discussed the prospects 
for the cosmological tests of such configuration. 

Finally, a conservative upper bound on type I seesaw right-handed neutrino masses, 
following from the requirement of electroweak vacuum stability, was derived, analysing in 
particular its dependence on the top mass. 
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Appendix A: Formulae for the RG running at NNLO 



1. The /3-functions 



Here we provide the expressions for the /3-functions up to 3-loops, see eq. (2). 
At 1-loop they are given by: 

rt\ 27 Q Q 

$ ] = -j-g{tf + -g'(t) 2 g(t) 2 - 9X(t)g(t) 2 + \g\tY - 36^(t) 4 + 4A(t) 2 - 3g'(t) 2 X(t) 

+l2h t {t) 2 X{t) , 

$ = §M*) 3 ~ l9(tfht(t) - 8g 3 (t) 2 h t (t) - l ^g'{tfh t {t) 
19 



& = , 
/?« = -7g 3 (t) 3 . 

At 2-loop they are: 

/3f = 80g 3 (t) 2 h t (t) 2 X(t) - 192g 3 (t)%(t) 4 + ^(i) 6 - ^V(f) W - yfc(t) W 

-^X(t)g(tf - *-fg'(tfg(t) 2 + Q3g'(t) 2 h t (t) 2 g(t) 2 + ^g'(t) 2 X(t)g(t) 2 - 3^(t) 4 A(t) 

+^^(t) 2 A(t)^(t) 2 - ^V(t) 6 + 180Mt) 6 - 16^(t) 2 ^(t) 4 - |A(t) 3 - y W 

69Q 8^ 

-24/^(t) 2 A(t) 2 + 6 {3g(t) 2 + g'{t) 2 ) A(t) 2 + ^g'(tfX(t) + ^^(t) 2 ^(t) 2 A(t) , 



if-' 



(2) 



23 



1 Q 

-108g 3 (t) 4 + 9g(t) 2 g 3 (t) 2 + -g\t) 2 g 3 (t) 2 + 36^(t) 2 ^ 3 (t) 2 



4 



9'(t) 2 g(tf 



4 9{t)A + W^ (t)4 - 12/lt(t)4 + + ht{tf 



^g{t) 2 + ^,'(t) 2 - 2A(t/ 



/3f = 12, 3 (t) 2 ^(t) 3 + {^g{tf + ^'(t) 2 - \h t {t) 2 ^ g{tf , 
# } = V(*) 8 + (^(t) 2 + ^ s'(t) 2 - ^(*) 2 ) ,'(t) 3 , 

$ = 93(t) 3 ( K l9(t) 2 -26g 3 (t) 2 + 1 jg\t) 2 -2h t (t) 2 ^ . 
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The leading terms in the 3-loop /3-functions of A and h t are [12]: 



266 1 ^99 

' + 32( 3 )g 3 (t)%(t) A + (-38 + 240( 3 ) g 3 (t) 2 h t (ty - (-5- + 36C 3 )^(t)' 



3 



(3) 



1 1 244 1 

- 48(Mt)%(t) 2 \(t) + -(895 - U96( s )g 3 (t) 2 h t (t)±\(t) + 
6 3 6 

1(11! _ 198C 3 )^(t) 6 A(t) + h-1224 + 1152( 3 )g 3 (t) 2 h t (t) 2 \(t) 2 + 
08 36 

1 1719 97 1 

-(— + 756C 3 )^(t) 4 A(t) 2 + Y4 h t {tf\{tf + ^(3588 + 2016C 3 )A(t) 4 
(-^ + 320C 3 )s 3 (t) 6 + - 114C s )|fe(f) W - ^0.(0 W 

+( if + f Cs)Mt)6 + ^(t) 2 ^(t) 2 A(t) + f ^(t) 4 A(t) + |/i,(0 2 A(0 2 - ^A(t) 3 

where C 3 = 1.20206... is the Riemann zeta function. 

The complete 3-loop /3-functions for the gauge couplings are [11]: 

= 3 -^9(t) 7 + 39,(t) 5 , 3 (t) 2 + 81*(t) W + | 5 (t)V(t) J - \g{tfg,{t) 2 g'{t) 2 
~9{tfg'{tY - 7 -§g(t)%(t) 2 - Utfg 3 {tfh t {tf - ™g(tm) 2 h t (t) 2 

^ = ^(*) VW S - 9(t) 2 9s(tmt? + 9% 3 (t)V(t0 3 + ^</(f) V(*) 5 - ™gs(t) 2 g'{tf 
388613 j 7 785 2 , 3 2 29 2 , . 3 2 2827 , , 2 



<?W - ^f)V(*w*r - - g3 (ty g '(tyh t (ty - — g \ t yh t (ty 



5184 J w 32" w " w w 3 ^v/-v/ *w 28g 
= X* (t) W + 21^(t)^ 3 W 5 + f^ 3 (t) 7 - ^(*) W<A*) 2 + y<fe(f) V(*) a 

-^*(*)V(t) 4 - f g(ty 93 (tyh t (ty - m 93 {ty ht {ty + ^ 3 (*)V(t) W 

+lhg 3 (t)%(ty. 



2. Higgs quartic coupling matching 

According to Sirlin and Zucchini [27], the 1-loop matching is given by 

«g»M = (f/iW + AM + ^) . (AD 
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where £ - 

A(aO 

/o(m) 



Z\z\ 
A(z) 



§■ and, introducing c 
3 



my 
m z ' 



2 

3c^ 

l>2 



log(0 - log (c 2 ) + 6 log 



m 



H 



^ log (c 2 ) + 12 log c 2 (c 2 ) + ^ 



-6 f 2( 2 
+2Z 



— ^ + 1 ) lof 

mi 



VTLr 



log 

3m 2 



+ 4c 2 Z 



9 /25 
~9~ 
15 



7T 

71 



(A2) 



- y (2c 2 + |^3) 



rrir 



4 log 



VTLr 



+ 2Z 



(A4) 



= 8 (2c 4 + 1) - 12c 4 log (c 2 ) - 12c 4 Z 



+ 6 2c 



4mf , 
+ 1 I lo f 



m, 



-6Z 



24mT / I m. 
+ I log 



m 4 ? \ \ /// 
2A{z) arctan (^j^y 
log )+1 



+ z 



- 2 



-A(z) 



i£z>\ 
if z<| 



(A5) 

(A6) 
(AT) 



We compute the QCD and the Yukawa contribution to X^(i^) following the expressions 
of [14] (multiplied them by a factor of 6 to compensate for the different definition of the 
quartic coupling). 



Appendix B: Seesaw contribution to the /3-functions 

Below the right handed neutrino mass scale, the running of the effective light Majorana 
neutrino mass is given by [49] 

dmAt) k (-3g 2 (t) 2 + 6h t (tf + ^f] m u (t) . (Bl) 



dt v 6 

For /i > M u , we have [50] 

^ = - Kit) (h v {tf + \h t {tf - \g'(tf - 9 -g(t)^ , (B2) 

together with 

6^ = -3/^(t) 4 + 2\(t)h u (t) 2 , 5$ = l -K{tf . (B3) 

2 
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